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I. INTRODUCTION 



The first result on the Bethe ansatz solutions of the asymmetric simple exclusion process 
(ASEP) is to be probably attributed to the short remark by Dhar [1]. The detailed calcu- 
lations were later given in a seminal paper of Gwa and Spohn [2], where the issue about 
the scaling of the next to the largest eigenvalues of the evolution operator was addressed. 
The subject has attracted a significant attention since then and a number of the exact re- 
sults has become available. Among the results obtained with the Bethe ansatz, there are 
the scaling function describing the crossover of the next to the largest eigenvalue from the 
Kardar-Parisi-Zhang to the Edwards- Wilkinson regime [3], the large deviation function for 
the distance travelled by a particle [4], [5]. These results describe the behaviour of the sys- 
tem in the large time limit. For the finite time evolution, we first of all should mention the 
Schiitz's determinant solution of the master equation for the ASEP on the infinite lattice 
[6]: the conditional probability of an arbitrary particle configuration, given an initial particle 
configuration, was obtained. This solution was later generalized to the ring geometry [7], to 
the ASEP with particle dependent velocities [8] and used to evaluate the distribution of the 
current of particles across the arbitrary bond [9] at the infinite lattice. 

All the above mentioned results concern the continuous time ASEP, where the jump 
of a particle is described by the Poisson process. The continuous time solution can be 
straightforwardly extended to the discrete time dynamics with random sequential update 
[5]. The Bethe ansatz solution for the ASEP with backward ordered update [10] was given 
in [11]. The history of studies of the ASEP with fully parallel update rises probably from the 
90s, when it was widely discussed as a simplest version of the Nagel-Schreckenberg traffic 
model with maximal velocity of cars: Vmax = 1 [12]. In that case the approximate mean- field 
method proposed to find the stationary measure of the traffic flow models [13] resulted in 
the expression, which becomes exact in the thermodynamic limit. Later it was realized that 
the mapping to the zero range process (ZRP) allows one to represent the stationary measure 
of the parallel ASEP at the ring in form of the product measure [14]. In the case of open 
boundary conditions the stationary measure was obtained in form of the matrix product 
[10], [15], [16]. At the same time, the description of the time evolution of the ASEP with 
parallel update remained challenging unsolved problem for a long time. Recently, it has 
been shown that the parallel ASEP is a particular case of the general stochastic process 
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with parallel update, which is solvable by the Bethe ansatz [17]. 

In present article we use the Bethe ansatz to obtain the solution of the master equation 
for the ASEP with parallel update on an infinite lattice. Given an initial configuration, we 
obtain the probability of an arbitrary configuration of particles. The solution is expressed 
in terms of the determinant of the matrix, which depends on the initial and final configura- 
tions of particles and time. Though the general idea of the solution is similar to that of [6] , 
the parallel dynamics leads to peculiarities, which make a generalization of previous results 
nontrivial. Below we give two independent derivations of the final result using two different 
treatments of the Bethe ansatz: the analytic representation based on the expansion of the 
solution over the Bcthc ansatz eigenbasis and the geometric one based on the direct enumer- 
ation of particle trajectories. While the former reveals the peculiarities of the cigcnspace of 
the evolution operator, the latter clarifies combinatorial aspects of the problem. Of course 
both methods lead to the same final result. 

The article is organized as follows. In the Section II, we formulate the model and sketch 
the results of the article. In the Sections III and IV we derive the solution using the analytic 
and geometric Bethe ansatz respectively. In the Section V we summarize the results and 
discuss perspectives. 

II. FORMULATION OF THE MODEL AND RESULTS. 

Let us consider P particles on ID infinite lattice. Below we imply that the particles move 
from left to right, which will be also referred to as a forward direction. At each step of 
the discrete time each particle may take one step forward with probability v or stay with 
probability (1 — v) provided that the next site is vacant. If the next site is occupied, the 
particle stays with probability 1. All sites are updated simultaneously at each time step. 
We define a configuration X of particles by a set of their coordinates X — {xi, ^2, . . . , xp}. 
Below we imply them to be ordered: 



The probability Pt{X) for the system to be in a configuration X at time t obeys the Markov 
equation 



— oo < Xi < X2 < ■ ■ ■ < Xp < oo 



(1) 




(2) 



{X'} 
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where T{X, X') is the probabihty of the transition from X' to X for one time step. The 
transition probabihties T{X, X') are defined by the above dynamical rules. 

The goal of the article is to show that the conditional probability, P{X, 0), for the 
system to be in a configuration X = {xi, X2, ■ ■ ■ , Xp} at time t, given it was in a configuration 
X' = {a;? ,X2,---, x^} at time 0, is the ratio of two determinants 

P(X tlx' 0) - detF(X,X°,t) 

P{X,t\X,0)- ^^^p^x^xQy (3) 

The matrix elements of the matrix F {X, Y, t) , which depend on two configurations X — 
{xi,X2i ■ ■ ■ 1 Xp} and Y = {7/1, 7/2, ... , yp} and on time t, are 



(F(X,y,i)).. = /(i-j>,-y,-,i), 



(4) 



where the function / (a, 6, t) is expressed in terms of the Gauss hypergeometric functions: 



f{a,h,t)^{l-vn-lfl 



( V \b i-t-a)b TP 
/ 



a, — t — a + 



V 



6+1 



' v-l 



(-Ml 2 



a — b, —t — a 
-6+1 



' v-l 



6 > 



b<0 



I 



(5) 



The notation (a)^ means the shifted factorial (a)^ = a(a + 1) • • • (a + n — 1). 



III. ANALYTIC APPROACH. 



A. Hilbert space and the vector of state. 

Let us introduce the Hilbert space supplied with the complete left and right bases con- 
sisting of vectors {X\ and \X) respectively, where X runs over all particle configurations, 
with inner product 

{X\X')^5{X,X'). (6) 
The state of the system at any time step can be associated with the state vector 

|P,) = J]P,(X)|X). (7) 

{X'} 
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In terms of the state vectors, the master equation takes a simple operator form 

|P,+i) = T|P,), (8) 
where the evolution operator T is defined as follows 

T= \X)TiX,X'){X'\. (9) 

The conditional probability P(X, t|X°,0) can be represented as the matrix element 

P(X,t|X°,0) = (X|T*|X°>. (10) 

To evaluate matrix elements, we construct the set of left eigenvectors \Bz) of the operator 
T 

T |P^) = A(Z) |P^) (11) 
and the adjoint set of right eigenvectors {Bz\ 

{Bz\T^A{Z){Bz\, (12) 

which, as will be seen below, are parametrized by the P-fold parameter Z. We will prove 
that the set forms a complete basis, i.e. provides the expansion of the identity operator 

J^{X\Bz){Bz\X')^{X\X'). (13) 
z 

As a result, the transition probability (10) we are looking for can be reduced to the evaluation 
of the sum 

{X\ T* |X°> = Yl <^|T*|5z> <^z|^°> = ^^1^^) (Bz\X') ■ (14) 

z z 

B. Bethe ansatz. 

The solution of the eigenproblem (11) corresponding to the master equation of the general 
integrable ASEP-ZRP model with parallel update was given in [17]. Here we sketch the 
solution for the particular case of the parallel ASEP defined above. 

For further purposes it is convenient to introduce a notion of cluster of particles. A 
sequence of n occupied sites situated one after another with two empty sites on the ends 
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without any other empty sites between them will be referred to as a cluster of the length n. 
A single isolated particle is the cluster of the unit length. During the evolution an isolated 
particle can approach and join a cluster from behind resulting in decrease of the number of 
clusters at the lattice by one. The opposite process is when a cluster splits up by a first 
particle detaching the cluster ahead creating an isolated particle and increasing the number 
of clusters by one. If a particle jumps between two clusters, which are separated by only 
one empty site, the number of clusters does not change. 

The explicit form of the transition probability T{X, X') is a product of fac- 
tors, each corresponding to a particular cluster of particles in the initial configu- 
ration X' . The value of these factors is either v or {'^ — v) depending on whether 
or not the first particle of a given cluster jumps during the transition from X' 
to X. 



T{X,X')^ Yl {1 - vf-"^' v""' (15) 

i=l 

Here mj = 0, 1 is the number of particles hopping from i-th cluster of X'. 

To proceed with the solution, we introduce a transformation [17] given by the diagonal 
in the configurational basis operator V 



\X) 'A1 



where the wight W{X) is the stationary measure of the process. For the fixed number of 
particles at the infinite lattice, it can be given in terms of the number of clusters He {X) in 
configuration X 

W (X) = (1 - . (16) 

Technically, it is more convenient to solve first the right eigenproblem 

To\B%) = K{Z)\B%) , (17) 

for the operator Tq defined as follows 

To = VTV^ 
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FIG. 1: The transformation of the transition probabiUties. 



rather than for the operator T, which have the same eigenvalue but different eigenvectors. 



Then the eigenvectors \Bz) of T can be related to eigenvectors of Tg 



\Bz)=V-'\B%). 



(18) 



Explicitly the matrix elements 



To(X,X') = ^^T(X,X'). 



(19) 



W{X) 

can be constructed in terms of the cluster structure of final configuration X according to 
the following rule: 

// a cluster of particles in X receives a particle from behind, it contributes 
the factor v to Tq{X,X') and (1 — v) otherwise. 



Nc{X) 

To{X,X')= H {l-vY'^'v'^ 

i=l 



(20) 



where fcj = 0, 1 is the number of particles approaching i-th cluster of X from behind. Indeed, 
the event when an isolated particle joins a cluster, the first particle of which stays, contributes 
f (1 — v) to T{X, X'). At the same time, as the two clusters merge into one, it follows from 
(16,19) that the additional factor of {1 — v)~^ should be incorporated into Tq{X, X'), so 
that the final contribution will be v (see Fig 1). If the first particle of a cluster of the initial 
configuration jumps forward, while no particle is received from behind, the probability of 
the jump V enters into T{X,X'). However, as one cluster in the initial configuration has 
split up into two clusters in the final one, the factor (1 — v) should be added in Tq{X,X'). 
The processes, where the both first and the last particles of the cluster stay, proceed without 



change of the number of clusters Mc {X), and hence their contributions to T{X, X') and to 
Tq{X, X') are equal as well as when the cluster of the initial configuration receives a particle 
from behind while its first particle detaches it at the front. The expression (20) can also be 
obtained in more formal way using the correspondence between the ASEP to the ZRP as it 
was done in [17]. As a result, (17) can be written in the form 

J^c(X) 

K{Z){..., x,x + l,...,x + n, ,...\B%) = E n ^'^(1-^)'"'' 

i-th cluster i'^*} 

(. . . , X - /Cj, X + 1, . . . , X + rij, . . . \B%), (21) 

where the index i in the product runs over all the clusters of particles of the configuration 
X. The number ki takes on the value of 1 or depending on whether or not the i-th cluster 
receives the particle from behind. The summation is performed over all possible sets {ki}. 
The i-th cluster consisting of Ui particles is written explicitly in the argument in formula 
(21). 

As usual, the strategy of the Bethe ansatz solution is the following. We note that when 
all the particles arc at least one empty site apart from each other, they behave like noninter- 
acting, i.e. each particle decides whether it is going to jump independently of the location of 
the others. In other words, if the coordinates of particles {xi, . . . , xp} of the configuration 
X satisfy the inequality {xi -|- 1) < Xj+i for every i, obeys the following equation 

A(z)(xi,...,xp|s°> = ^n^^'(i-^)^-'^^ 

{ki} i=i 

x{xi-ki,...,xp-kp\B''z). (22) 

Here the product is taken over all P particles, and the summation is over all ki, which take 
on the values of 1 and depending on whether or not the i-th particle decided to jump. The 
necessary trick is to find a way for this equation to be satisfied in the whole set of particle 
configurations including the domains where the equality {xi -|- 1) = Xi+i takes place. If we 
consider formally the equation (22) for such a configuration, the terms like {. . . ,x,x, . . . IB^) 
will appear. They are beyond the allowed domain (1) where the sets of particle coordinates 
are defined. Therefore they should be redefined and cancelled so that the difference with 
(21) disappears. Particularly, one needs to make sure that the result of the one cluster 
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summation in (21) coincides with the corresponding term of the free equation, 

V . . , X — 1, X + 1, . . . , X + n, . . . \Bz) + {1 — v) . . , X, X + 1, . . . , X + n, . . . \Bz) 

n 

= J2Il'"''*i'^-^f''''{---^^-ko,x+l-k2,...,x + n-kn,...\B''z). (23) 

{ki} i=0 

This relation can be satisfied with the only constraint on B'^ {xi, . . . ,xp) 

V [{..., x-l,x,... |5|> -{...,x-l,x+l... |5|>] 
+ {l-v)[{...,x,x,... > -{...,x,x + l,... \Bl)] = 0, (24) 

which can be checked by a direct calculation for the size of the cluster n = 2, and then 
proven by induction for larger n (see [17]). 

After that one can look for (xi, . . . , xp) in the form of the Bethe ansatz, 

(xi, . . . , xp|i?i> = 5^ ^.,....p^.T • • • (25) 

parametrized by the parameter Z, which is a set of P complex numbers Z — {zi, . . . , zp}. 
Here the summation is performed over all permutations a of the numbers 1, . . . ,P. The 
substitution of (25) to the free equation (22) yields the expression for the eigenvalues in 
terms of parameters {zi, . . . , zp} 

p 

A{Z) = il-vfl[{l + Xz,), (26) 

i=l 

where 

A = — ^. (27) 
1-v ^ ^ 

The substitution of the Bethe ansatz into the constraint (24) results in the relation be- 
tween the amplitudes A^-^^^^-p, which differ from each other only by two neighboring indices 
permuted 

^^-S{z,,z,), (28) 

where 

To obtain an explicit form of the amplitudes A^^, consider a particular permutation 
((7i, . . . , (7p). Let us associate the factor ^j^7^ to an elementary transposition (. . . . .) — > 
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{. . . j,i . . .). Apparently the power of will then coincide with the position of i in the 
permutation. Thus, we have the following correspondence 




(30) 



From (28), we have 



and 



= (-i)"'""!! (-^^1^)"". (32) 

where V ({c}) is the parity of a permutation {a}. The formulas (18,25,32) together with 
the definition of the weight (16) is all what we need to write the right eigenvector \Bz) as 
a set of its coordinates {X\Bz) in the configurational basis {|X)}. 

To construct the adjoint set of the left eigenvectors {Bz\ one has to solve (12) or equiv- 
alently the right eigenproblem for the transposed matrix T* 

T* \Bz) = A{Z) \Bz) . (33) 

To this end, we note that in the original process corresponding to the matrix T, the first 
particle of a cluster contributes the factor v to the transition probability T{X, X') if it 
jumps forward and {1 — v) if it stays. In the adjoint process, which corresponds to the 
transposed matrix T* the particles move in the opposite direction. The last particle of a 
cluster contributes the factor v or {1 — v) to a matrix element if it jumps approaching the 
cluster or stays in it respectively. One can see that the action of the operator T* is similar, 
up to the inversion of the coordinates, to the one of the operator Tq, (19). In other words, 
the matrices Tq and T* are related by the identity 

T* = ITqI (34) 

where I is the operator of the inversion of coordinates, I \xi^ . . . , Xp) = \—Xp, . . . , —Xi). It is 
then straightforward to show that the solution of the equation (33), corresponding to the 
eigenvalue (26) is given by the Bethe ansatz 

{Bz\X) = ^J2^,^-^^z:i...zZ (35) 
10 



with the amphtudes obeying the relation 

^^-S{z„z,). (36) 
^...ji... 

where the 5"- matrix is the same as in the relation for A^, (29), with the permuted arguments. 
As a result we obtain 

■^ai...ap — l/^<Ti. ..(Tp 

The constant factor (P!)^ is introduced for further convenience. Finally {Bz\X) and 
{X\Bz) can be represented in the determinant form 

{X\Bz) = W^(X)detB, (37) 
{Bz\X) = -^detB (38) 

where the matrix elements Bij and Bij of the matrices B and B respectively are given by 

""-'/^'^-{-^r^r- (39) 

C. Completeness of the eigenbasis and the transition probabilities. 

To proceed, we should define the range of the summation over Z in (13,14). For the 
infinite lattice, the spectrum of Z can be taken continuous, and the summation over Z 
should be replaced by the P— fold contour integral 

The contour of the integration over each Zi is closed around the points z = and 2; = 1, so 
that the point z — — 1/A stays outside. Fig. (2). This is a generalization of the solution for 
the continuous time ASEP. In the limit A — > the point z — —1/A goes to infinity, and the 
situation becomes similar to the one studied in [6] . Below, we prove the completeness of the 
eigenbasis, which follows from (13). As soon as it is established, one can directly substitute 
(25,26,32) into the expression for the probability (14) and obtain 

p 

{x\ T* |x°> = w{x) J2 i-^r Uf(' - - <' ^) (41) 

{a} i=l 
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FIG. 2: Contour of integration over the parameters zi, . . . , zm- 

where (—1)°^ is the sign of the permutation a and the function /(a, b, t) is defined as follows 

/KM)^(l-.)'/^(l + A.,'(-i±^)".- (42) 

We should note that initially there should be a double sum over permutations a and a' as 
the terms {X\Bz) and {Bz\X) contains one sum each. However, there is no any dependence 
on a' in the product of the factors f{a[ — a^' ,x^'_ — , , t) over all is. As a result we have 
P\ equal terms under the sum over a', which just eliminate the prefactor 1/ (P!) in (35) and 
the only one sum (41) over permutations remains. 
Consider two configurations 

X = {xi, < X2,< . . . ,< xp} ,Y = {yi, < 7/2, <•••,< yp} • (43) 

The explicit form of (13) is given by the particular case of formulas (41,42) taken at t = 0, 
= Y. Consider a term of the sum corresponding to a particular permutation a, which is 
the product of P integrals of the form. 

/(.-.....-,,..0)^/^(-i±^)'""'.--,. (44) 

Apparently, the summand to be nonzero, the integrals for all is should be nonzero simulta- 
neously. The following statement can be proved. 

Proposition 1 The expression (44) yields a nonzero result for all i = 1, . . . , P only when 
X = Y and only for those permutations, which permute the particles inside their clusters. 

It means that the relations 

Xi = yi, (45) 
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Xi 



(7,;. 



(46) 



hold for any i = 1, . . . , P . 

We should remark that in the case of the continuous time ASEP, which corresponds 
to the particular limit of our model with A ^ 0, [6], the contribution to the continuous 
time version of (13) comes from the only summand, which corresponds to the identical 
permutation ai — i. In the Appendix A, we prove the above proposition for the general case 
under consideration, while the continuous time case follows from the proof. 

Given the configuration X = Y with any pair of coordinates Xi, x^i satisfying (45,46), we 
can calculate the integral (44), which is equal to (— A)*^°"' for i > (Ji and 1 for i < The 
summation over all permutations of particles within the clusters leads to the determinant 
of P X P matrix 



j dZ{X\Bz) {Bz\X) = W{X) detFo, 
which has a block diagonal form 



(47) 



S(ni) 











S(n2) 



S(n3) 



(48) 



Here, S(n) is the n x n matrix 



S(n) 



1 -A (-A)^ (-A)=^ 

1 1 -A (-A)' 

1 1 1 -A 
111 1 



(49) 



which corresponds to a cluster consisting of n particles. It can be transformed to the 
triangular form by subtracting neighboring columns, which allows the calculation of its 
determinant: 

detS(n) = (1 + A)"-' = (1-^;)-"+' (50) 

The product of such terms over all clusters yields just the value inverse to the weight W{X) 
defined above: 

det Fo = (1 - t;)^=W-^ = 1/W{X) (51) 
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The substitution to the formula (47) completes the proof of the relation (13). We should 
note that the matrix (49) has a lower triangular form in the limit A ^ 0, with all the units 
on the main diagonal, as should be in the continuous time limit. 

The integral (42) can be evaluated in terms of the hypergeometric functions. The sum 
in the r.h.s. of the formula (41) can be rewiten in the determinant form 

{X\ T* \X') = W(X) det [fix, - ^r?, i - j, t)] , (52) 

while as follows from (51) the factor W{X) is the inverse of the determinant 

W{X) = 1/ det [fix, -xj,i- j, 0\^^,_M (53) 

The resulting expression is given in the Section 2, (3-5). 

Once the answer has been obtained, one can directly reexamine it. First, we note that 
the correct initial conditions 

detF(X.XO.O)^ 

detF(X,X,0) ^ ' ^ ^ ^ 

follow from the Proposition 1. We should note that unhke the continuous time case, the 
contribution to the determinants comes not only from the elements of the main diagonal of 
the matrix F{X, X, 0), but from the blocks connecting the particles within the same clusters. 
Then the equality 

detF(X,X",t+l) ^ ^ detF(X^X^t) 

detF(X,X,0) ^ ' MetF(X',X',0) ^ ^ 

can be proved in spirit of the proof given by Schiitz [6]. We first use just established 
fact (53) that the determinants in the denominators in both sides are the inverse stationary 
weights 1/W{X) and 1/W{X'). Due to (19), brought to r.h.s., they transform the transition 
probability T(X,X') to To(X,X') 

detF(X,X°,t+ 1) ^^TQ{X,X')deiY{X',X^,t). (56) 

{X'} 

Now one should use the properties of the matrix elements f{xi — x'jji — j, t) to show that the 
equality (56) holds. The solution being the hnear combination of the eigenvectors, all the 
properties of its constituents follow from the corresponding properties of the Bethe ansatz. 
Particularly the constraint (24) imply that the function f{a,b,t) obeys the relation 

= [A/ (gi, ni-l,t) + f (qi, Hi, t)] [f {q2, n2, t) - f (q2, n2 + 1, t)\ 

- [A/ (g2 - 1, 712 - 1, + / (^2 - 1, ^2, t)] [f (gi + 1, m, t)-f {qi + l,m + 1, t)] (57) 
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where ni,n2 are arbitrary integers, qi — 1 — P, . . . , P — 2 and q2 — 2 — P, . . . , P — 1. As 
above this relation ensures that in the physical domain the master equation can be written 
in form of the equation for noniter acting particles. The latter is satisfied if every individual 
function /(xj — a;°, i — j, t) for any pair i,j satisfies the equation for a free single particle 

f{q, n,t+l)^ vpfiq, n-l,t) + {l-v) f{q, n, t). (58) 

Both relations can be checked by substitution of the explicit form of the function /(a, b, t) 
and follow from the contiguous relations for hypergeometric functions. Thus, the expression 
obtained is indeed the sohition of the master equation. 

We should note that if one would use the resulting determinant form as an ansatz with an 
unknown function f{a,b,t), the relations (57,58) can be considered as difference equations, 
which are to be solved to find this function without the above Bethe ansatz analysis. Such 
a representation can be considered as an alternative treatment of the Bethe ansatz. Of 
course it leads to the same final result. A problem, however, occurs on the way. To solve 
the difference equations one should impose the initial conditions, e.g. /(I — P, n, 0), which 
reproduce the correct value of the resulting determinant. Up to our knowledge there is no a 
regular procedure of doing this. In fact we encountered the equivalent problem in the above 
Bethe ansatz analysis, which is how to choose the measure of integration over Z. Luckily, 
the simplest choice made yields the correct result. We, however, can avoid considering this 
problem at all by using the constructive approach to the probability derivation based of the 
direct enumeration of particle trajectories. This approach is presented in the next section. 

IV. GEOMETRIC APPROACH. 

In this section, we give an alternative solution of the ASEP with the parallel update based 
on a combinatorial analysis of particle trajectories. We introduce a geometrical treatment 
considering sequentially three models of interacting particles with increasing complexity of 
interaction. The simplest one is the free-fermion model of non-intersecting trajectories, 
known in physical literature as the problem of "vicious walkers" [18] and in mathematical 
literature as the Gessel-Viennot theorem [19]. Then, we consider the ASEP with the back- 
ward sequential update to demonstrate how the Bethe ansatz appears with a complication 
of interaction. Finally, we apply the developed approach to our main problem of the ASEP 
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with the parallel update. 



A. Combinatorial ansatz for vicious walkers 

Consider P particles labelled {1, 2, . . . , P} hopping in one direction on an infinite one- 
dimensional lattice. The discrete space-time dynamics of their motion can be described 
as a set of continuous broken trajectories on a triangle lattice A, which is obtained from 
the square lattice by adding a diagonal bond between the upper left corner and the lower 
right corner of each elementary square. Let {x, t) be integer space-time coordinates of a 
particle on A, where the vertical time axis is directed down and the horizontal space axis is 
directed right. A trajectory of each particle is a sequence of connected vertical and diagonal 
bonds of A. The diagonal bonds correspond to jumps of particles to their right for a unit 
time-step with probability v. The vertical bonds correspond to stays of particles at fixed 
sites during the unit time-step with probability 1 — v. The initial positions of particles are 
Xi<X2<---< x%. Let our particles be found at sites xi < X2 < ■ ■ ■ < xphy the moment of 
time t. The problem of vicious walkers is to find the conditional probability for the particles 
to reach the positions X = {xi,X2, ■ ■ ■ ,xp);xi < X2 < ... < xp from the initial positions 
x^,X2, ■ ■ ■ , iCp) for time t so that no pairs of trajectories are intersected during time 
t. In other words, one assumes that every site can be occupied by at most one particle at 
every moment of discrete time and if this rule is violated at a moment f < t, the process 
stops. 

We start with consideration of the one-particle motion on the infinite lattice. Let %{x\x^) 
be a set of one-particle trajectories, which are starting at 0) and finishing at (x, t). Each 

trajectory q G Tt{x\x^) is realized with probability v^^'^"{l — vy--^+^" , Therefore, the total 
probability for the particle to reach x from x^ for time t is 

P{x,t\x^,0) =v''-'\l-vy-^+^"\\Tt{x\x'')\\ =Fo{x-x^\t). (59) 

where 

Fo{x\t) = (%^{l-vy-^ (60) 



yXy 

is the weight of one-particle trajectories in the set Tt{x\x^). For the case of P particles, the 
set §p of all possible free trajectories is a direct product of one-particle sets of trajectories 
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reaching {X,t) from {X°,0): 

Sp = Ttixilx'l) Tt{x2\xl) ® . . . ® %{xp\x^p) (61) 

The set Sp contains non-intersecting and intersecting trajectories. The latter cases are 
forbidden by the single-occupation rule and should be subtracted in evaluations of probability 

P{X,t\X^,0). 

Consider first the case P — 2. We denote the set of forbidden elements by U12, U12 C §2 
where indices 1, 2 emphasize that the order of final points of intersecting trajectories coincides 
with the initial one. To cancel the contribution of forbidden elements, we introduce an 
auxiliary set of pairs of trajectories 

A2i^Tt(x2\xl)®%{xi\xl), (62) 

where the final coordinates are permuted with respect to those in §2- It is easy to see that 
all elements of the set A21 are pairs of crossing trajectories. Each intersecting pair q & S2 
has a first collision point {xc, tc), i.e. the space-time point, where the trajectories meet for a 
first time. Consider the set S2(xc, tc) of intersecting trajectories with the fixed first collision 
point {xc,tc). The sets S2(xe,tc), (a^c^c) G A break the set U12 into subsets parameterized 
by coordinates {xc,tc). For all {x'^jt'c) 7^ {x",t") we have 

S2«,i;)n§2«,iD = (63) 

and 

y S2(x„g=Ui2 (64) 

For every set E)2{xc,tc), there exists a uniquely defined set A2i{xc,tc) C A21 obtained 
from E>2{xc,tc) by permutation of tails of all trajectories beginning at the first collision 
point. The sets E)2{xc,tc) and A2i{xc,tc) are geometrically equivalent, as there is one-to- 
one correspondence between their elements. Therefore ||S2(a;c, tc)|| = ||A2(a;c, tc)||- Like 
the sets E>2{xc,tc), the sets A2i{xc,tc), {xc,tc) G A break the set A21 into subsets: for all 
«,^c) 7^ «>^c) we have 

A2,{x'cX)^^2i{x:X)-^ (65) 

and 

U A2i(a;e,g = A21. (66) 
{xc,tc)eA 
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It follows from (64,65,66) that the whole sets U12 and A21 are geometrically equivalent 
and IIU12II = ||-'^2i||- Given {X,t) and (X°,0), all elements of S2 have the same weight 

2 

i=l 



Then, we have for the probabihty 

P{X,t\X',0) = g X (||§2|| - IIA21II) = (67) 

= (Fo(a;i - x[\t)Fo{x2 - a;°|t) - Fo{x2 - x[\t)Fo{xi - xl\t)) = 
= detM, 

where the elements of the 2x2 matrix M are 

M,,, =Fo(a;,-a;°|t), i,j = l,2. (68) 

This result may be easily generahzed for the P-particle system of vicious walkers [18]. 
The set §p contains the subset of intersecting trajectories Ui2...p. We introduce P! — 1 
auxiliary sets 

= Tt{x^^ \xl) Tt{x^, \xl) . . . ® Tt{x^^ 14) (69) 

where a — {ai, . . . , ap} is any permutation of numbers 1, 2, P beside the identical one. 
Each element of the set 

Ui2...pU{U^-> (^0) 

containing trajectories intersecting at given point {xc,tc), has a unique geometrically iden- 
tical counterpart with a pair of permuted indices, e.g. 12. ..i...j...P and 12...j...i...P. Intro- 
ducing the sign of permutation (—1)'^, we can ascribe opposite signs to two counterparts. 
Then, we obtain the known result [18]: 

P{X, t\X', 0) = g X I^IISpII + ^(-1)-|| A,|| j = (71) 
= detM, 

where 

Mij = Fo{xj-x^,\t), t,j = 1,2,. ..P (72) 
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B. Combinatorial treatment of Bethe ansatz 



The vicious walkers are locally interacting particles. In this section, we consider the sys- 
tem of particles with the interaction of a non-zero range, which takes place in the case of 
the totally asymmetric exclusion process with the backward sequential update [15]. Con- 
sider again P particles hopping in one direction on an infinite one-dimensional lattice. The 
interaction between particles can be defined by the following rules: 

1. Trajectories of particles do not intersect (every site can be occupied by at most one 
particle) . 

2. A particle stays at its own site with probability 1 if the target site is occupied by 
another particle during the step of discrete time. 

Like the previous case, the trajectory of each particle is a sequence of vertical and diagonal 
bonds on the lattice A. Each diagonal bond has weight v and the vertical bond weight 
1 — V. In view of more complicated interaction, it is convenient to decompose the set of all 
free trajectories into more detailed subsets. For each vertical bond [{x,t'), {x,t' + 1)], the 
trajectory of i-th particle passing this bond can be decomposed into two parts 0) — > 
{x,t')\l\{x,t' + 1) {xi,t)] and [(x°,0) {x,t')\ - v\{x,t' + 1) {xi,t)] where the value 
between vertical bars means the weight of the bond [{x, t'), {x, t' + l)]. These new trajectories 
are geometrically equivalent, but the first trajectory passes the selected vertical bond with 
weight 1 and second one with weight —v. We make this decomposition for each vertical 
bond of each trajectory and denote the whole set of decomposed trajectories of i-th particle 
by Tt{xi\x'^) using the same notation as for the set of one-particle trajectories in the previous 
section. 

The weight of the trajectory q G Tt{xi\x^) is a product of the weights of its bonds. 



The weight of the set of trajectories is a sum of weights of its elements. In accordance with 
the above definition (59), the weight of the whole set %{xi\x^) is 



t 



ljL{q) — y^fj,{k—th. bond of q) 

k=l 



(73) 



fi{Tt{xi\x^J) = Fo{x^ - x^i\t) 



(74) 



The set of free P-particle trajectories §p. 



§P = Tt{xi\x^i) ^ %{x2\xl) . . . (g) %{xp\x%) 



(75) 
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has the weight 



p 




(76) 



qeSp i=l 



As in the case of vicious walkers, the set Sp contains a subset of unallowed elements,Ui2...p, 
which should be excluded. By the rules 1 and 2, an element of Sp is unallowed, if there is 
at least one pair of intersecting trajectories, or if there are two neighboring vertical bonds 
from which the left one has weight —v. To cancel unallowed elements we start as above with 
the case P — 2 and construct an auxiliary set of trajectories A21. We will see that, due to 
nonlocality, the auxiliary set has more complicated structure 

00 1 

A21 = U U - h) ® Tt{xi\xl - k2). (77) 

ki=0 k2=0 

Beginning the construction, we notice that every unallowed pair has a first collision point 
{xc,tc) G A, where the particles come for the first time to neighboring sites at a moment 
tc < t. The first trajectory reaches the site {xc,tc) from (x5,0) and second one reaches the 
site (xc + 1, tc) from (xg, 0) and then it makes the vertical step to the site {xc + 1, tc + !)• 

The first trajectory has, just after the collision, a diagonal bond [{xc, tc), {xc + l,tc + 1)] 
with weight v (refered to as a collision of the first type) or vertical bond [{xc, tc), {xc, tc + 1)] 
with weight —v (refered to as the coUision of the second type). Notice, that trajectories, 
which have the vertical bond [{xc, tc), {xc, tc + I)] with weight 1 are allowed. Thus, we have 
two types of unallowed elements of U12 for fixed {xc,tc). Denote these subsets by Y{xc,tc) 
for the first type and W{xc,tc) for the second one. For all {Xc,t'J ^ {.^'ci't'c) we have 



For each space-time point {xc, tc) G A, we construct the sequence of pairs of trajectories 
•^(^1,^2) and A^(A;i,A;2), (^1 = 0,1,2,... and /c2 = 0, 1, arguments Xc,tc are omitted). 



w{x'c, t'c) n w{x'^, tl) = 0, 



and 




(78) 
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obtained from V and W by permutation of tails in each pair of trajectories and shifting 
the initial parts in negative direction. 

Explicit form of transformations V =^ A^{ki,k2) and W ^ A^(fci,/c2) are given in 
Appendix B. Using these relations, we can express the probability 

P(X,t|X0,0) = M§2\ U12) (79) 
via the weight of all unallowed configurations 

{xc,tc) k=0 

Our next aim is to bring this expression to a determinant form similar to (72). Consider 
operator a^, which acts on the set of the free trajectories of i-th particle and gives the set of 
free trajectories with the shifted origin by one lattice space to the left: 

aiTt{xj\x^i) = Tt{xj\x^i - 1) (81) 

The operator representation allows us to write (80) in compact form: 

Pi{ U12) = (j^Jt{x2K) (8) Tt{xi\xl)^ , (82) 

where the denominator is defined by its expansion 

^ 00 

^ = E«?- (83) 

^ fc=0 

and the summation of operators means the joining of the sets. 

The crucial property of the operator expression in (82) is its factorization with respect 
to indices 1 and 2. Introducing the functions 

00 

Fi(x2 - xl\t) = /X ((1 - ai)-^%(x2\xl)) = Yl ^0(^2 - X? + k\t), (84) 

fe=0 

and 

1 

F_i{xi-xl\t)=fi{{l-a2)Tt{xi\x'2))=J2Fo{xi-xl + k\t) (85) 

fe=0 

we obtain the two-particle probability P{X, t\X^, 0) = //( S2) — A*( U12) in the form 

P{X, 0) = Fo{xi - x^i\t)Fo{x2 - 4|i) - Fi{x2 - x°|i)F_i(xi - 4|i) (86) 
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or 



P{X,t\X°,0) ^det M, 



(87) 



where 




(88) 



To generalize the determinant formula (87) to the general P-particle case, we should take 
into consideration several additional properties of intersecting trajectories. First, we may 
organize the procedure of exclusion of unallowed elements of U12...P in an ordered way. 
Starting with the top line of the lattice A, we examine all free trajectories of the set §p, 
row by row, until we meet the first collision point where unallowed trajectories are cancelled 
with elements of the auxiliary set A^-. 

If the number of particles P > 2, the elementary squares associated with the collision 
point {xc,tc), {Xc,tc){xc + l,tc){xc,tc + l){xc + l,tc + 1), may occur several times in one 
horizontal strip of A. If squares filled by interacting trajectories are separated one from 
another by a gap of empty sites, the above arguments can be applied to each pair of inter- 
acting trajectories independently. The crucial case for the Bethe ansatz is a situation when 
the elementary squares are nearest neighbors. The specific property of the totally ASEP is 
that, in each pair of interacting trajectories, the right trajectory remains free and interacts 
with the next trajectory independently on its left neighbors. Therefore, we can analyse the 
interaction between particles considering successively elementary squares in each row from 
left to right starting from an arbitrary empty square until all unwanted trajectories will be 
removed. 

After t' steps from the top to bottom, one obtains the set of path configurations q E Sp 
which are allowed in the first f rows, and the set of auxiliary configurations yet not involved 
into cancellation procedures. Remembering that all elements of A^- for all cr ^ 1 have end 
points permuted with respect to the original order 12... P, we conclude that each element of 
Ag. contains at least one collision point. Therefore, all elements of the auxiliary set will be 
cancelled after t' — t steps with unallowed elements of U12...P. 

The described way of exclusion of unallowed configurations implies a successive construc- 
tion of the auxiliary set A^. We notice that each two intersecting trajectories in A21 are 
non-equivalent: one of them belongs to the particle which overtakes another and can be 
called "active". On the contrary, the second particle is "passive". In the case P > 2, one 
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trajectory can overtake m others, and we call it "m-active". Similarly, the "m-passive" 
trajectories appear. 

Assume, that the trajectory of a given particle has m active intersections. It means that 
it participates m times in the cancellation procedure and its starting point is shifted m times 
to arbitrary distances in the negative direction. As a result, the auxiliary set associated with 
the free trajectory between points x° and Xj becomes 

Ttixjlx"",) (89) 



(1 - fli) 

Similarly, for trajectories having m passive intersections we get 

(l-dir%{x,\x^) (90) 

The weights of sets of m-active and m-passive trajectories are given by functions intro- 
duced in [6]: 

-— -— 7;(x,|x°) ) = Yl ~ ^o{xj - x° + k\t), (91) 



\ 



and 



F_m{xj - x'^lt) = ((1 - d,r%{xj\x'^)) = E I I (-'^)'Fo{x, - x° + k\t), (92) 

fc=0 \ ^ / 

Activity m of each trajectory is defined uniquely by the permutation of numbers a, so 

we have for the weight of auxiliary set A^- 

p 

M A,) = J]F,,_,(a:,,-x°|i) (93) 

i=l 

This product is a term of expansion of the determinant detM with matrix elements for 

Mij = Fi_j{xi - x°|i) z,j = 1, 2, P (94) 

It follows from (75) that the term corresponding to the identical permutation cr = 1 is 

p 

M§p)=J]Fo(a;,-x°|i) (95) 

1=1 

Collecting the contributions from the set Sp and all auxiliary sets A^., we obtain the 
determinant formula [6] 

P(X,i|X°,0) =detM, (96) 

which is vahd for all P > 1. 
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C. Combinatorial solution for parallel update 

The TASEP with parallel update is characterized by the rules 

1. Trajectories of particles do not intersect. 

2. A particle stays at its own site with probabihty 1 if the target site is occupied by 
another particle just before the step of discrete time. 

The first rule coincides with that for the the backward sequential update, however the 
second rules are different. The combinatorial solution for parallel update accumulates all 
main steps of two preceding sections. We start with the construction of different types of 
unallowed elements of U12 at the fixed collision point {xc,tc)- Two of them, Y{xc,tc) and 
W{xc, tc) are the same as for the sequential update. Two new sets X(xc, tc) and Y{xc, tc) 
are 

X(a;c, tc)'- the first trajectory has, just after the collision, the diagonal bond [{Xc, tc), {Xc + 

1, tc+ 1)] with weight v and the second one has also the diagonal bond [{xc + l,tc), {xc + 

2, tc + 1)] with weight v; 

Y{xc, tc)' the first trajectory has, after the collision, the vertical bond [{xc, tc), {xc, tc+1)] 
with weight —v and the second one has the diagonal bond [{xc + 1, tc), {xc + 2,tc+ 1)] with 
weight V. 

Like (78) we have 

(J ( Y{xc, tc) U W{xc, tc) U X{xc, tc) U Y{xc, tc)) = U12 (97) 

(a;c,tc)GA 

Similarly, we introduce along with the auxiliary sets Ay{ki,k2) and A^[ki,k2) new sets 
Ax{ki, k2) and Ay{ki, k2) obtained from X and Y by permutation of tails in each pair of 
trajectories and shifting the initial parts by integers ki and /c2- Explicitly, the transforma- 
tions X ^ Ax{ki, /C2) and Y =^ Ay{ki, k2) are given in Appendix B. 

The whole auxiliary set A21 has, instead of (Bl), a more complicated structure. 
To reveal it, we show schematically in Fig. (3) the positions of sets V, W, X, Y and 
At,(0,0), A^(0,0), Ax{0,0), A^(0,0). It is easy to see that all these sets are character- 
ized by reference points from where trajectories reach the end points xi and X2 after the 
collision. For the set V, the point xi is reached from Xc + 1, the point X2 from Xc + 1, for 
the set W, Xi is reached from Xc, X2 from Xc + 1; for the set X, xi is reached from Xc + 1, 
X2 from Xc + 2; and for the set Y, xi is reached from Xc, X2 from Xc + 2. The coordinates 
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FIG. 3: The reference points and the weights of the sets V, W, X,Y (full circles) and 
A^(0,0),A^(0,0),A^.(0,0),A^(0,0) (full squares) 

of the reference points of V, W, X, Y are shown in Fig. (3) as full circles. 

For the auxiliary sets Ay{0, 0), Au,(0, 0), Ax{0, 0), A^t,(0, 0), the coordinates of reference 
points can be obtained by permutation of xi and X2- The positions of auxiliary sets are 
shown in Fig. (3) as full squares. 

The weights of sets in the "scattering zone" between tc and + 1 are: v{l — v) for 
V,-f(l -v) for W, for X, -v^ for Y,v{l-v) for A^(0,0), -v{l - v) for A^(0,0), 
w^) for A^.(0,0) and -t;2) for Ay(0,0). 

Our aim is to combine the sets Ay{ki,k2), Aw{ki, k2), Ax{ki,k2), Ay{ki,k2) in such a 
way that all sets of unallowed trajectories are cancelled with their auxiliary counterparts. 
Using the shift operators di, we consider sequentially four operator forms: 



1 



1 + Xd^^ ' 



(98) 



A, 



1 



1 ~ di 
is = 1 + \d^^ 



(99) 



(100) 



A, 



1-02 



;ioi) 



where A = v/{l — v). 

Each shifted set contains trajectories where an initial part from t = Otot = tc + l 
is geometrically equivalent to that of non-shifted trajectories. It means that the reference 
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FIG. 4: The transformation of the reference points due to action of operators Ai {a),A2 (b),and 
AsA^ (c) 

points are also shifted under action of dj. As a result, the action of operator Ai leads to the 
diagram Fig. (4a), the action of operator A2 leads to the diagram Fig. (4b), and the action 
of operator gives the final diagram Fig. (4c) where positions of auxiliary sets coincide 

with those of unallowed sets V, W, X, Y together with their weights. 

Thus, we may conclude that the whole set of auxiliary trajectories can be obtained by 
action of operator y4iyi2^3^4) so we have 

(l-a2)(l + Aar'; 



/i(Ul2) = /i 



(102) 



(1 - di) (1 + Adg 

which should be compared with the analogous expression (82) for the sequential update. 

The operator expression (102) is also factorizable with respect to indices 1 and 2. It 
allows us to introduce functions which are generalization of Fm and F_m given by (91,92): 



and 



or, in the explicit form: 



F^Nt) = J2Y1 fj}^ ^^ u (t^ ^0 AT + kt), 

^-^ ^-^ A; + n !n! m — n ! 

n=0 k=-n V ' ^ V I 



and 



n=0 k=—n 



[k + n)\nl{m — n)\ 



(103) 



(104) 



(105) 



(106) 
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The remainder of derivation coincides with that for the sequential update and we can write 
for P{X, 0) the determinant formula with the matrix elements 

Mij = Fi_j(xi - x'jlt) i,j = 1, 2, P (107) 

The determinant formula for P(X, i|X°,0) is valid for all initial states and for all final 
states X conditioned by inequalities Xj+i — > 1 for alH = 1, 2, P — 1. If Xj+i — Xi — 1 
for some 1 < i < -P — 1, the determinant formula should be corrected. To find the correction, 
consider the final positions of particles Xi = ^2 = x + 1 at time t in the two-particle case 
P = 2. For the last time step, the final positions are simultaneously the reference points of 
unallowed and auxiliary sets of trajectories. However, these sets cannot be superposed as in 
the case of bulk scattering because the range of auxiliary set exceeds the range of unallowed 
set (the points {x,x + 2) and {x + l,x + 2) in Fig. (4) are absent). Exphcit calculations for the 
transitions from (x — 1,^ — 1), {x + l,t — l) to {x,t), {x + l,t) and from {x,t — l), {x + l,t — l) 
to {x,t), {x + l,t) show that the correct probability differs from the determinant result by 
the factor (1 — v). Using the fact that the scattering of each pair of trajectories can be 
considered independently, we write the final result as 

P(X,t|X°,0) = (1 -i;)"detM, (108) 

where n is the number of pairs of neighboring particles in the final position. The equivalence 
of the results (3-5) and (108) is verified directly by expanding the hypergeometric functions 
in (5). 

To make above arguments more transparent, consider a transformation X X^ of a 
realization of the process interchanging the initial and final configurations of particles. As 
a result, we obtain the same realization overturned (Fig.5a,b). 

The initial positions in Fig. (5a) and final positions in Fig. (5b) obey the conditions Xj+i — 
Xi>l for all i — 1,2, P— 1. The final positions in Fig. (5a) and initial positions in Fig. (5b) 
contain pairs of neighboring coordinates. A remarkable property of the path configurations 
is that the difference in numbers of bonds having weight 1 (yellow bonds in Fig. (5a) and red 
bonds in Fig. (5b)) equals to the difference in numbers of pairs in final and initial positions. 
Therefore, we can apply the pure determinant formula to the case when the final position is 
free of neighboring pairs and then calculate the opposite case multiplying the determinant 
by (1 -t;)" as in (108). 
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(a) (b) 



FIG. 5: (a) The initial state contains no pairs of neighboring sites; the final state contains one pair, 
(b) The reversed dynamics. The yellow and red bonds have weight 1. The difference between their 
numbers equals to the number of pairs in the final state in (a) 

V. DISCUSSION AND CONCLUSION. 



To summarize, we considered the discrete time asymmetric exclusion process with parallel 
update at the infinite 1-D lattice. We have presented two alternative derivations for the 
probability of a particle configuration at arbitrary time, given an initial configuration. In the 
first approach we constructed the transformation of the evolution operator, which allows one 
to find its left and right eigenvectors in the Bethe ansatz form. We proved the completeness 
of the eigenbasis, which allows calculation of the matrix elements of an arbitrary power of 
evolution operator. The second approach is based on the consideration of the ensembles of 
trajectories of noninteracting particles. Taken with alternating signs their statistical weights 
are chosen in such a way that they add up to reproduce correctly the dynamics of interacting 
particles. 

The result for the system of P particles is expressed in the form of the determinant of 
P X P matrix, whose elements are given in terms of the Gauss hypergeometric functions 
2-Fi(a, b, c; x). To relate this result to the previous ones obtained for the ASEP, we note that 
in the continuous time limit t — oo, A — > 0, tA = r the hypergeometric function 2-Fi(a, b, c; x) 
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degenerates into the confluent hypergeometric function yielding 



f{a,b,r) = {-ire-- 




iFi (a-6, -6 + 1,t) 6< 



iFi {a,b + l,T) 



b>0 



(109) 



which is equivalent to the results obtained by Schiitz, [6]. The relation of Schiitz's functions 
to the confluent hypergeometric functions was already noticed before [20]. 

We should also remark that there is a direct relation of our results to the determinant 
solution for the TASEP with backward ordered update [11] referred also to as a fragmentation 
model [9]. The determinant formula for that process can be obtained from our one by 
using (1 — v)"'f{a,b,t — a) instead of f{a,b,t) in the formula for the matrix elements (4). 
This connection between conditional probabilities seems to be similar to the correspondence 
noted in [9] on the level of the current distributions. Given such a close relation between 
the solutions of the master equations, an establishing of the direct mapping between the two 
original processes would be interesting. 

The obtained results open rich perspectives for derivation of more complicated correlation 
functions for the density and current of particles. Many of such quantities have been obtained 
recently due to close connection of the subject with the statistics of random permutations and 
ensembles of random matrices. Particularly, using the mapping of the ASEP with parallel 
update to the growing Young diagrams Johansson obtained the probability distribution 
of the integrated particle current through an arbitrary bond starting from the half filled 
lattice [21]. Properly rescaled this distribution converges to the distribution of the largest 
eigenvalue of the random matrix of the Gaussian unitary ensemble [22] . The limiting case of 
the Johansson's formula for the continuous time ASEP as well as its analogue for the ASEP 
with the backward update were reproduced using the Bethe ansatz in [9] . A derivation of the 
Johansson's formula starting from our determinant expression can be proposed as a closest 
goal for further investigation. These results were generalized for other initial conditions, 
which turned out to be related with different ensembles of random matrices [23], [20]. In 
addition more intricate quantities were considered like multipoint joint distribution of the 
integrated current [24], [25] or correlation functions of site occupation numbers [26]. Most of 
results were obtained either for continuous time ASEP or for the ASEP with the backward 
update. One can hope that the results of present paper can serve as a starting point for the 
further studies of the case of parallel update. 
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Another direction of research is the generahzation of the results for the ring geometry. 
In the present article, the proof of the completeness of the Bethe ansatz basis for evolution 
operator of the parallel ASEP has been given for the case of the infinite chain. Similar proof 

of the completeness and orthonormality for the case of the ring would be an important issue 
not only for the stochastic processes community but in the general context of the Bethe 
ansatz solvable models [27] . 
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APPENDIX A: PROOF OF THE PROPOSITION ??. 

The integral (44) can be evaluated using the Cauchy theorem. To this end, one can 
expand the expression in the parentheses into the Laurent series. Practically one can first 
expand (1 + Az)'"^* and {1 — 1/ zY^~^ separately and then evaluate their product. Then the 
coefficient of the term of power [xi — yo-j) will be the result. The choice of the integration 
contour means, that the expansion of (1 + \zy~'^^ will contain only nonnegative powers of \z 
and the expansion of (1 — l/z)'^*"* will contain only nonnegative powers oil/ z. In addition 
we should take into account that the number of terms of the expansion will be finite when 
the exponent is positive. Thus, when 

i = (Ji, (Al) 
the integral to be nonzero we should require 

Xi = Ui. (A2) 
In the other cases, the following inequalities should hold. 

i-Ui > Xi-y^i, for i > Ui (A3) 
(7i-i > Vai- Xi, for i < Ui. (A4) 
Moreover the domain (43), where the coordinates of particles are defined, implies 

i-(Ti< Xi-Xai, i - (Ti < yi - y^^, for i > ai (A5) 
CTi-i < x^i- Xi, (Ti - i < y^i - yi, for i < ai. (A6) 
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Comparing the expressions (A3,A4) with (A5,A6) respectively, we obtain 

Vai > x^i, Hi > Xi for i > Ui (A7) 
Va < XcTi, Vi < Xi for i < Ui. (A8) 

Any permutation can be decomposed into the set of disjoint orbits [28], i.e. the sequences 
of elements which trade places with one another under a given permutation. Consider an 
orbit ii —> ^2 ?fe ii, such that 

C^ii = «2,(7i2 = ^3, • • • ,(7j^ = ii (A9) 

There exists r such that ir is the largest in a given orbit. Then from > ir+i = and the 
second inequality of (A7) we have 

Vir > Xir, (AlO) 
while from i^-i < ir = o'i^_-^ and the first inequality of (AS) it follows that 

Vir < (All) 

These two are consistent only if 

Vir = Xir- (A12) 

The latter equahty together with the inequahty (A3) and the second inequality of (A5) yield 

Vir - Vir+l >ir- V+l > Vir " Vir+l^ (Al3) 

which means 

ir - ir+i = Vir - Uir+i- (A14) 
In the same way using the inequalities (A3,A5) we obtain 

ir ir—1 — Xi^ -^v+i ■ (A15) 

We should note that the equality like 

Xi — Xj = i — j (A16) 

means that there are as many sites between the sites Xi and Xj as many particles. In other 
words i-th and j-th particles belong to the same cluster. Apparently, the coordinates of any 
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pair of particles between them also satisfy similar relation. One can see that the particles at 
sites belong to the same cluster of the configuration Y, while the particles at sites 

Xi^,Xi^_^-^ belong to the same cluster of the configuration X. Let ir+i be the largest of the 
numbers v+ij V-i- The cluster of particles spreads from Xi^ backward to at least Xi^_^_-^ in 
both configurations X and Y. Therefore we conclude that 

Xi,+^=yir+i- (A17) 

Then we look for ir+2- If V+2 is larger than ir+i, then it is between ir and i^+i and hence 
belongs to the same cluster in both X and Y and wc should look for ir+s- If V+2 is less than 
ir+i, using the same arguments as for ir, v+i we have 

Vir+i - yir+2 = *r-+2 - ir+i, (A18) 

i.e. we extend the cluster in Y backward up to ir+2-th particle. Then we again compare ir+2 
with ir-i, which currently limits the cluster in X from behind. If i^-i is between ir+i and 
ir+2 it necessarily belongs to the same cluster also in Y. Then the cluster can be extended 
up to ir-2, etc. Thus, step by step we extend the cluster backward in Y going forward along 
the orbit and backward in X going backward along the orbit until we come to the same 
element in both directions, which means that the orbit is exhausted. Finally we conclude 
that all the particles in a given orbit should belong to the same cluster and their coordinates 
in the configurations X and Y coincide. 



APPENDIX B: EXPLICIT FORM OF AUXILIARY SETS 

The transformation V =^ Ay{ki,k2)is 

[(x°,0) ^ (xc,tc) {xc + l,tc + l) (xi,t)]0[(xl,O) (x^ + l,tc) (xc + l,tc + l) 
(X2, t)] [(x? - ki, 0) ^ (xc - ki, Q ^ (xe + 1 - /ci, + 1) ^ {X2, t)\ ® [{xl - /ca, 0) 

(Xc + 1 - /C2, tc) {Xc + l- k2, tc + 1) ^ {Xl,t)], 

and the transformation W Aw{ki, k2) is 

[(X?, 0) ^ {Xc, Q {Xc, tc+1) ^ {Xi, t)] (g) [{xl 0) ^ {Xc + 1, tc) ^ {Xc + 1, tc + 1) 

{X2,t)] =^ [(x? - ki,0) {xc - ki,tc) {xc - ki,tc + 1) ^ {X2,t)] (g) [{xl - k2,0) 

{Xc + 1- k2, tc) {Xc + 1- k2, tc + 1) ^ {Xl, t)]. 



32 



It follows from the definitions that 



oo 



U U U iMkl,k2)U A^{ki,k2))^ A21. 



(Bl) 



(xc,tc)eA ki=0 k2=0 



The weights of introduced sets obey the following relations 



/x( A„(0,0)) 



(B2) 



//( A,(0,1)) 
Ai( A^(/c,0)) 



li{A^{k + l,0)), k 



0,1,2... 



(B3) 
(B4) 



II { A^{k,l)) 



li{A^{k+l,l)), k 



0,1,2... 



(B5) 



Transformation X =^ Aj.(A;i, /C2) is 

[{Xl, 0) ^ (Xe, tc) {Xc + 1, + 1) (Xi, t)] ® [{xl 0) ^ {Xc + 1, tc) ^ + 2, + 1) ^ 
(Xa, t)] ^ [K - /Ci, 0) ^ (Xc - /Cl, tc) ^ (Xe + 1 - ki, tc + l) ^ {X2, t)] ® [(x^ - /Cs, 0) ^ 

{xc+1- k2, tc) ^ {Xc + 2- k2, tc+1) ^ {xi, t)] 
and the transformation Y =^ Ay{ki, /C2) is 

[(x?, 0) ^ (xc, tc) (xc, tc + 1) ^ (xi, t)] (8) [{xl 0) ^ (Xc + 1, tc) ^ (Xc + 2,tc+l) ^ 

{X2,t)] [(x? - /ci,0) ^ (xc - ki,tc) ^ (xc - A;i,tc + 1) (a;2,t)] ® [(x^ - k2,0) 
(Xc + 1 - k2, tc) ^ (Xc + 2 - k2, tc+1) ^ (xi, t)]. 
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